function dydx = test(x,y)
global e1;
dydx = [y(2)
    50^2*(-e1-1+x^2/2)*y(1)];